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Abstract 



Static observers remain on Killing- vector world lines and measure the rest-mass+kinetic 

g; ^^^^^^^ 

space and time. The total mechanical energy is the total flux through a spacelike cut at 
one time. The difference between the total mass-energy and the total mechanical energy 

o 

is the total gravitational energy, which we prove to be negative for certain classes of sys- 



tems. For spherical systems, Misner, Thorne and Wheeler define the total gravitational 
energy this way. 

To obtain the gravitational energy density analogous to that of electromagnetism we 

X 

first use Einstein's equations with integrations by parts to remove second order deriva- 
tives. Next we apply a conformal transformation to reexpress the scalar 3-curvature of 
the 3-space. The resulting density is non-local. 

We repeat the argument for mechanical energies as measured by stationary observers 
moving orthogonally to constant time slices like the "zero angular momentum" observers 
of Bardeen who exist even within ergospheres. 
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1 Introduction 



In classical physics, energy has different forms which are additive. Energy is conserved and 
one can assess how much is transferred from one form to another. In Einstein's theory 
of gravitation different forms of energy in the matter are mixed in inseparable ways with 
gravitational binding energy. Gravitational energy is diffuse, partly mixed up with various 
forms of matter energy, partly stored in the gravitational field itself. Alas, the matter is the 
source of gravity so that it is impossible to disentangle gravitational energy from other forms 
of energy with or without solving the field equations. This is a great loss that adds to the loss 
of "gravitational force" and "gravitational energy density" with which it shares a common 
origin: the principle of equivalence. 

Misner, Thorne and Wheeler 1 [32] gave an expression for the total gravitational binding 
energy of a spherical system. They also suggest a possible definition for the gravitational 
energy density in this particular case. Katz J3] realised that their definition of gravita- 
tional energy could be extended to all stationary systems and gave explicit expressions for 
gravitational energy density of stationary systems whose spaces were conformally flat. 

Here we give an expression for the gravitational energy density of any stationary space 
that is asymptotically flat. Our expression is non-local in that it involves the conformal 
transformation that removes the scalar 3-curvature of the stationary 3-space and it is invariant 
under transformations of spatial coordinates. However in specially chosen coordinates (such 
as isotropic ones when they exist) it may be expressed locally in terms of derivatives of the 
metric components. 

MTW's starting point is the difference between the Schwarzschild mass found asympoti- 
cally at large r and the total mechanical energy which is the sum of the rest mass + internal 
+ kinetic energies of the matter. This difference is the gravitational energy. Kinetic en- 
ergy has to be measured relative to some standard of rest and even when that standard 
is agreed it will not be conserved; as fluid climbs out of gravitational wells kinetic energy 
decreases etc. However as fluid moves from one place to another there will be a flux of 
mechanical energy and in relativity we consider flux vectors describing how the mechanical 
energy moves from one region of spacetime to another. In any stationary spacetime we may 
better known as MTW.... 
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use the time-like Killing vector 2 to define 'static' observers whose velocities w 11 = 
(£ is the magnitude of are unit vectors along Such observers assess the mechanical 
energy of a mass mo moving with 4- velocity u^ 1 to be vtlqc 2 / 'y/l — v 2 /c 2 = itiqcPu^w 11 , v is 
its velocity with respect to the local observer. Likewise a fluid of dust of proper rest-mass 
density p (in its own rest frame) will have a rest-mass energy density in any frame moving 
with 4-velocity n p of pc^u^n^ = pc 2 /\/l — v 2 /c 2 because the dust particles seem closer to- 
gether due to Lorentz contraction of their proper volume. Such factors occur in the following 
expression for rest mass. For a fluid of dust the rest-mass flux vector pu^ is conserved. The 
total rest-mass-energy is given by 

M c 2 = [ pc^u^^/^gdT.a where cE M = -re^oadx" A dx p A dx a , (1.1) 
js 3! 

and the 3-surface £ is any space-like hypersurface that spans all space. To get the mechanical 
energy flux vector associated with this fluid as assessed by static observers we must weight the 
rest-mass flux vector with the mechanical energy per unit rest-mass 1 / a/1 — v 2 /c 2 = u^w^ . 
We obtain 

pu^u v w v ^gdY,^ = T?w v 'v^cEm- ( L2 ) 
The total mechanical energy on any cut through spacetime will be 

E M = J^w v ^dT l ^. (1.3) 

We have motivated expression (|1.3|) above by using a fluid made up of dust because the 
argument is easiest to follow in that special case; however the final expression in term of Tjf 
is not at all confined to a dust fluid nor even to a fluid at all. The same expression would 
still hold if part of the total T]f were the Maxwell stress energy-momentum tensor for the 
electromagnetic field. Likewise it would still hold for a plasma with an anisotropic pressure 
tensor or any other physical field. 

Whereas formula 1)1 -lj) holds for dust on any cut through spacetime, formula (|1.3|) will 
give answers that depend on the cut chosen unless Tjfw 1 ' is a conserved vector which is not 
generally so. That said, there are many important special cases for which it is conserved. 
Since 

D M (T>") = Dp = -T>"0 M log& (1.4) 

2 Indices A, fj,, v, p, ■ • ■ = 0, 1, 2, 3; indices k,l,m,n - ■ ■ — 1,2, 3, the metric has signature H and g is 

its determinant. Covariant derivatives are indicated by a D, partial derivatives by a d. And for once neither 
G nor c are set equal to 1 and a; = ct. Finally the permutation symbol in 4 dimensions is e^ vpa with £0123 = 1 
and in 3 dimensions by e feim with £123 = 1. 
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whenever the flux vector T^w v lies in the equipotential surfaces of constant £, it is conserved. 
In stationary spacetimes with purely toroidal motions in (t, x 1 , x 2 , <p) coordinates, Tjfw u is 
only non zero with [i in the t or (p components whereas <9^1og£ is only non zero with jj, in 
the other components so in all such cases T{fw u is conserved. When that occurs, we may 
evaluate Em over any cut through spacetime and all will give the same answer. Generally 
we may always evaluate Em but should not expect to get an answer independent of the cut 
chosen. Thus the mechanical energy seen by static observers can depend on the cut through 
spacetime over which it is evaluated. In practice there is normally a good time coordinate 
such as Boyer-Lindquist in Kerr and cuts are chosen to be at constant time. 
Following [23] the total gravitational energy is 



where Mc 2 is the total energy. In this paper we convert the gravitational energy into "field 
energy" like Maxwell's (E 2 + B 2 )/8ir in electromagnetism. To do this we first reexpress the 
T/f in Em by using Einstein's equations. Retaining the scalar 3-curvature we remove other 
second derivatives via integrations by parts. The scalar 3-curvature is reexpressed employing 
a conformal transformation. Thus the total gravitational energy Eq is found in terms of 
gravitational field variables. 

So far our whole argument has been developed for kinetic energy as seen by the static 
observers following the Killing vector field but such observers do not exist inside ergospheres. 
Furthermore Bardeen [2] has emphasised that the static observers have angular-momentum in 
spaces with spin. He has introduced zero-angular-momentum observers (ZAMOs) who move 
along the normals to the hypersurfaces of constant time. Although such observers move 
around in azimuth they may be preferred to our static observers because the latter give 
zero motion to a static body which has angular-momentum (backwards). Furthermore these 
observers exist even within the ergosphere. Thus, t-hypersurface orthogonal observers move 
with a different 4-velocity field so they assess a different kinetic energy than the static 
observers. According to them the mechanical energy is not Em but 



When the mechanical energy is assessed over a hypersurface of constant time, ()1.6|) seems 
far more natural than our expression (jl.3j) : for dust the u^w^ = — v 2 /c 2 factor is the 

same as u M n M because the normal n M to the hypersurface E is w^. Writing Tjf = pu^Uy in the 
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Eq = Mc 2 - E M 






above expression thus yields {p/y/1 — v 2 / c 2 ){l / 1 — v 2 / c 2 ) integrated over all space. The 
first '7' factor comes from the Lorentz contraction and the second from the mechanical energy 
per unit rest-mass. The factors are now the same - they are not normally the same for static 
observers because for rotating systems their worldlines do not have orthogonal hypersurfaces. 
Thus there are strong arguments for preferring the hypersurface orthogonal ZAMOs to the 
static observers and for considering 



as the gravitational energy rather than Eq. In axial symmetry the ZAMO mechanical en- 
ergy splits into two terms l E — as in classical mechanics. We discuss the angular- 
momentum density of gravitational fields elsewhere. 

In what follows we derive the gravitational energy densities for both sets of observers. 
These densities are different and integrate to different totals. We shall naturally assume that 
the dominant energy conditions hold flU so that, in particular, 



A proper evaluation of gravitational energy may be useful in numerical modelling of 
relativistic stars [37], relativistic ellipsoidal configurations JSj an d thick spherical shells 
among others. Limits of stability are characterized by the relative binding energy (M — 
Mb) /Mb, Mb is the "baryonic mass", which is taken as a relativistic generalization of the 
T/|W| ratio used in astronomy |6; • However, the more relativistic the configurations, the less 
representative that ratio is because (M — Mb)c 2 is a mixture of all forms of energy and in 
some models that "binding" energy even changes 3 sign! The T/| W| ratio is better generalized 
as Ek/Eq where Ek is the bulk kinetic energy as measured by local observers. 

There are other properties of gravitational energy which will not be dealt with here but 
are interesting. For instance, it is not difficult to show that if gravitational energy for static 
observers is extremal (5Eq = 0), spacetime must be flat. It might be possible to show, along 
the lines used by Brill and Deser [7J, that Eq, near flatness, is a local maximum. It may be 
true but is as yet unclear that gravitational energy as defined here is always negative. We 
show that Eq is indeed negative for static and some stationary systems. 

3 A typical example is the binding energy of Einstein's spherical shell of self-bound particles |13| (see also 
flip^l. A recent discussion on the sign of the binding energy in static perfect fluids can be found, for instance, 
in Karkowski and Malec |26|. 



E G = Mc 2 - E M 



(1.7) 



T ilv w^w v > and T^w^Wy > 0. 




5 



2 Stationary spacetimes 



Here is a short summary of some basic formulas for stationary spacetimes in which we also 
fix notations. A good recent summary of present knowledge about stationary spacetimes is 
given in Beig and Schmidt's paper 4 [3]. 

(i) Coordinates adapted to stationary spacetimes, and static observers 

Stationary spacetimes have a timelike Killing field We consider first static observers: 
they have velocity components in the £^ direction: 

w fl = j where f = yt^, (2.1) 

and are at rest with respect to fixed observers at infinity. 

In the projection formalism, the metric is decomposed as follows: 

ds 2 = (w x dx x Y + (g^u - w^w u ) dx^dx v . (2.2) 

In practical calculations the metric is most often written in coordinates in which it takes the 
following form HE], ^ , jHS]: 

ds 2 = g^dx^dx" = f(dx° - A k dx k ) 2 - lkl dx k dx l \ (2.3) 

in these coordinates, 

{e} = {1,0,0,0} and {^} = {/- 1 / 2 , 0,0,0}. (2.4) 

f,A k , and 7fc/ are functions of x k only. plays a special role: indices of tensors in 7-space 5 
will be displaced with that metric and its inverse. Thus, for instance, we shall write: 

A k = 1 kl A l and A 2 =A k A k = -f kl A k Ai>0. (2.5) 

The metric components g^ and their inverse g^ v are: 

floo = / , 9ok = -fAk , 9kl = ~1ki + fAkAi, (2.6) 
g°° = r l ~A 2 , g ol = -A l , g kl = -j kl . (2.7) 

4 See also |38|: the original formulation of the projection formalism is due to Geroch [16). 
5 Also called the "quotient space" obtained by quotienting spacetime by the action of the stationary isometry 
3 . It represents the collection of the orbits of the Killing vectors 

The 7fci give a measure of proper lengths and / of proper times. See a good discussion in |28| . See also |39|. 
Indices of and are always displaced with g^v not ju or later with ju. Thus in l|3.2()^ = gkvT hh ' . 
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The determinants g of g^ v and 7 of 7^ are so related: 

V=g = f 1/2 V^- ( 2 -8) 

The metric (j2.3j) and the components w p , Q2.4|) . keep the same form under a change of 
coordinates: 

x '° = x + r°(x fc ) and (2.9) 

= T fc (x ; ). (2.10) 

Ak and 7^; transform like tensors for x' k = T k (x l ). The interesting transformations are 
however those of x° that leave / and 7/^ invariant but A' k = Ak + 9fcT°. They are associated 
with different x° = const hypersurfaces which we need to fix if we want a definite gravitational 
energy density. The usual gauge condition is 6 

d k (VlA k )=0. (2.11) 

The dominant terms of the asymptotic solutions for isolated sources are given in MTW [32] 
on page 456. In properly chosen asymptotically Minkowski orthogonal coordinates at infinity 

/-►I , Ai -> -2 - , 7w->-%(H ) where r = J x 2 + y 2 + z 2 -> 00, 

(2.12) 

and 

n fc = - , £(n fc ) 2 = l > > i k = ^- ( 2 - 13 ) 

Mc 2 is the total energy of spacetime and J = {Jk} = {— J k } its angular momentum vector. 
Einstein's equations, 

1 &7lG 

RfMu = k(T^ u - - g^T) with k = (2.14) 

are better written in terms of projected components that are gauge invariant: 

R^vf , -w p w p ){g U(T -w u w a ) and R^W - w»w p )w v . (2.15) 

In these expressions in our coordinates, Ak must necessarily appear in a curl. Let represent 
covariant derivatives in 7-space; we shall write 

7 fc V; = V k and V fc V fc = V 2 . (2.16) 



6 The gauge condition 112.1111 reduces the freedom of translations to such t° that satisfy the conditon 
V 2 t° = [see 12.161 for the definition of V 2 ]. If we ask for r = const at infinity this fixes x° up to a constant. 



The gauge invariant components of the Ricci tensor, invariant for time translations (|2.9f) and 
covariant for spatial coordinate transformations (|2.1U|) . are equivalent to |28| 7 : 

#00 = \ V 2 /-^V fc /V fc / + / 2 V^V^], (2.17) 

R ki = g ^ g iu R = _-L V V/ + -^V k fV l f + 2f lrnr M k A m] V [l A n] + TZ kl , (2.18) 

2/ 4f z 

R k = g ku R u = —VifV^ - fVi(V^ k A\ (2.19) 
and the 4-scalar curvature 

R = f^Roo ~ lkiR kl - (2.20) 

In (|2.18|) . 1Z kl are the components of the Ricci tensor of the 7-space, not of the x° = const 
hypersurface. The following linear combination of (|2.17|) and (|2.18l) is of particular interest: 

- lkiR kl - r'Roo = -lkiK kl - = -11- 3/V^V^j; (2.21) 

7Z is the scalar curvature of the 7-space. The left hand side of (|2.21|) can be rewritten using 
Einstein's equations (|2.14|) : 

- lkiR kl - f^Roo = -2kT ° + 2A k R k . (2.22) 

Therefore, (|2~2Tj) and fl2~22]) imply 

- n - 3/V [ U' ] V [fc> An - 2A k R k = -2kT$. (2.23) 

The last two terms on the left hand side can be recast in another form using ()2.19|) : 

- 3fV^V [k A l} - 2A k R k = -V [fc (/^)V[ fc ^] - V fc {2fAM k A 1 ^) . (2.24) 

Inserting this new expression into (|2~25i then transferring -V fc (2fAiV [k A 1 ^ from left to 
right we get 

- K - V [fe (M q )V^ = -2kT ° + V fc (2/i ; V^l) . (2.25) 

Before going further, it is perhaps interesting to rewrite Einstein's equations in terms of 
"gravoelectric" and "gravomagnetic" vector fields E k and B m living in 7-space in the spirit of 
HH1, EDI, ESI and [10]. The gravoelectric potential 

^ = log V7, (2-26) 

7 In the 1959 edition of Landau and Lifshitz there is apparently a printing mistake in the third equation on 
page 301. What in their notations is written +-|/ Q7 //3 7 should be — ■|/ a7 //5 7 - 
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and the gravomagnetic vector potential is A k . £ k is denned by: 



S k = -d k ^p and -f kl £ k £i = £ 2 > 0. (2.27) 

B m is divergenceless 8 : 

(d k Ai-d l A k ) = r lklm B rn or B rn = rf nkl d k Ai with r, klm = ^e klm ; (2.28) 
in particular, 

V [k A^ k A l] = ^ VklmV kln B m B n = \ B m B m = \b 2 > 0. (2.29) 
In terms of £'s and B : s the components of the Ricci tensor take the following forms: 

f- l R 00 = -V-£ + £ 2 + ^fB 2 , (2.30) 
R ki = +V k £ i _ £ k £ i + 1 j (^ki B 2 _ B k B i^ + n ki^ ( 2 31 ) 

f- l R k o = ^r l klm (-V l B m + 3£ l B m ) = ^{-VxB + 3£xB} k . (2.32) 

Also, 

R = -2V • £ + 2£ 2 - - fB 2 - K, (2.33) 
and equality (|2.25|) takes this form: 

-K- f(^B 2 +B-Ax£^j = -2kT ° + V • (J A x B) . (2.34) 

Dots and cross products are defined in 7-space. 

(ii) Metric components adapted to observers with velocities orthogonal to constant time slices 

The coordinates are the same as in (i) but observers are different. They have different 
4- velocity components; we denote them by instead of w^. is associated with the 
1 + 3 decomposition familiar in the Hamiltonian formulation of Einstein's equations 9 . By 
definition, coordinates exist in which 

= if 1/2 , 0>0,0}. (2.35) 
The metric (|2.2j) with replacing now takes the following form: 

ds 2 = f{dx ) 2 - %i (dx k - W k dx°) (dx l - W l dx°) , (2.36) 

8 In US] and B's are replaced by H = f /2 B. This H is called B in |55| . 

9 The original reference is Dirac |12|. See also pQ. Modern versions are given in various books like |32| and 
frT|. The components W k in 12.371 correspond to the shift vector N k and / 1/2 to the lapse N. 
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in which W are the components of the local coordinate velocities in units of c: 

111 fl rpk 

W k = =^ = ^-r; also set W 2 = % t W k W l . (2.37) 



The metric components and their inverses are: 

5oo = f-W 2 , gol = W l = ^ lk W k , g k i = -% h (2.38) 
g 00 = J' 1 , g ol = f- l W l , g kl = -^ kl + f~ l W k W l . (2.39) 

The determinant g is related to the determinant 7 of 7^; like this: 

y/=9 = f 1/2 V*f- (2-40) 

The metric (j2.36|) and the velocity components (|2.35jl keep the same form under a change of 
coordinates (|2.9|) and ()2.10j) . Some of Einstein's equations are now better written in terms 
of Einstein's tensor components G^: 

The simplest forms belong to the projected components of G^ v , see (|2.15j) . with replaced 
by w 11 . In our coordinates, EQ 10 : 

2/G 00 = K + K 2 - K kl K kl , (2.42) 

Rh = n kl + v (fe A) - £k£i + 2f 1 / 2 K^ k v [l) w m] - r 1/2 £ m w m K M 

+ V m (r 1 / 2 W m K kl ), (2.43) 
~f! 2 Gl = JV 2 RI = -VAK{-5{K), (2.44) 



and 

R = 2(fG 00 -j kl R kl ). (2.45) 

In these formulas TZ k i is the Ricci tensor of the spatial metric %i; TZ = ^f kl TZ k i the corre- 
sponding scalar curvature while K k \ (the second fundamental form of x° = const) is: 

K k i = r 1/2 V (fc iy ; K = K k = ~ 1 kl K kl while S k = -d k log f 1 ' 2 . (2.46) 

Vfc is a 7-covariant derivative. The time coordinate must be fixed if we wish to integrate 
over a definite spacelike hypersurface. It is fixed up to a constant by a gauge condition like 



The reference contains the clearest and most detailed calculations. We rewrote the formulas in our 
notations. 
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1)2.11(1 . The most common choice in the Hamiltonian formulation and in numerical relativity 
are maximal spacelike hypersurfaces K = which imply that observers are expansion free, 
D^w^ = 0. In stationary spacetimes the condition reduces to 

V k W k = equivalent to K = 0, (2.47) 

as follows form ([2.46 1) . Notice that the condition is identically satisfied in axial symmetry 
[see (|4~Td|) below]. 

(Hi) Relations 

The metrics in the form ()2.3)) - ()2.7l) and (|2.36j) - (|2.39[) are the same metrics in the same co- 
ordinates, adapted to different decompositions associated with different families of observers. 
The /, W and jki are related as follows to the /, Ak and jki'- 

f , fA k 

and reciprocally, 

f = f-W 2 , A k = -f~ 1 W k , iki = lki + fAkAt. (2.49) 
Notice that A k = j kl Ai, not f l Ai. We now turn to the calculation of gravitational energy. 

3 Gravitational energy 

(i) As calculated by static observers 

If w x are the velocity components of static observers, the mechanical energy they evaluate 
on x° = const is, as pointed out in 1)1.3)1 . 

E M = j T^w v ^gdY, tl = f T$dV where dV = ^fd 3 x = (1 - fA 2 )~ 1/2 dV. (3.1) 

The right hand side is written in our special coordinates. dV is the proper volume element 
in 7-space and not of the x° = const hypersurface which is dV . These two coincide in static 
spacetimes. 

Now consider ()2.34l) . Divide by 2k, multiply by ^7 and integrate over the whole space. 

The boundary conditions imply that the integral 11 of j^dk (yfyfA x B) is zero. Notice that 

1 We use abundantly Stokes theorem |41| . or Gauss theorem I28|. or Gauss-Ostrogradsky theorem |19| or 
Green-Ostrogradsky theorem |17| or Green's theorem )27| or the divergence theorem Hadamard 12(11 gives 
it no name and regards it as a generalization of J f'(x)dx = f(b) — f(a). It is useful to remember that the 
divergence theorem applies to divergences of continuous and single valued functions. 
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if the gauge condition (|2.11l) was not fixed, j^dk (y/jfA x B) would not be fixed and the 
remaining integrand in ()3.2j) would also be gauge dependent. It follows now that the gravi- 
tational energy, see Q1.5|) . 

E G = Mc 2 - E M = Mc 2 - [ [K + /(- B 2 + B ■ A x £)]dV. (3.2) 

2jK J 2 

We shall next eliminate the Mc 2 from ()3.2|) and obtain an integral that depends entirely on 
field components and their derivatives. 12 To this effect, we make a conformal transformation 
of the 7-metric, say, 

7w = e2x 7fc/, (3-3) 



and set, by analogy with ([2.27)1 . 

Fk = -dkX; also 7 W ^ = JP^ = JF 2 . (3.4) 
The conformal scalar curvature 1Z* and 1Z are so related: 

e 2x n * = n + AV k F k - 2T 2 . (3.5) 
We then define x by the condition that TZ* = 0, that is 

U = -4S7 k T k + 2F 2 , (3.6) 



or 



V 2 e^ 2 - hl^ 2 = 0. (3.7) 

The 7-metric becomes spherical at large distance and all spherical metrics are conformally 
flat. Therefore we shall ask the factor e _2x to behave like the conformal factor of 7^/ near 
infinity as in (|2.12l) : 

_ 2v 2m m 
e x —> 1 H or x ~^ f° r r ~^ 00 • (3-8) 

Equation ()3,7j) has been considered by Cantor and Brill who found a necessary and 
sufficient condition for an asymptotically flat metric to be conformally equivalent to another 
asymptotically flat metric with zero scalar curvature. They showed that the condition for 
the existence and uniqueness of a solution of (|3,7j) is that 

j Vq-VqdV >-- / q 2 TZdV (3.9) 



12 Notice that on a general hypersurface in a general spacetime one may specify the lapse function J 1 / 2 
arbitrarily by arbitrary time reparametrizations x'° = T°(x a ). However, in stationary spacetimes there are 
privileged coordinates in which the metric has the form (12.311 and the allowed transformations of coordinates 
are those given by 12.911 and 12.1UH . Our gauge condition 12.1111 fixes the spacelike hypersurface. 
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for all smooth functions q ^ of compact support. Notice that (|2.23|) with Rq = kTq can 
also be written as 

K+pB 2 = 2Kf- 1 T 00 . (3.10) 

This shows that in non-static spacetimes outside matter 1Z = — ^fB 2 < and conditions 
()3.9|) are not automatically satisfied. 13 

If a solution of ()3.7|) exists, we may replace TZ in (|3.2|) by its equivalent ()3.6|) and notice 
that because of (|3.8|) 

— / Vi{A^)dV = - [ T l dS h = - (~) Anr 2 = -Mc 2 . (3.11) 

2K J K Jr^oo K \ r z J 

As a consequence of (|3.11|) . Mc 2 disappears from (|3.2|) and we obtain for Eq: 

Eg = -- [ {T 2 + -J[B 2 + 2B-(Ax £)}}dV = f e G dV; (3.12) 
n Jx° 4 J x o 

we used (|3.1|) to reexpress dV in terms of the proper volume element dV on the slice x° 
= const. With the gauge condition (|2.11|) . the gravitational energy density ec on the fixed 
x° = const hypersurfaces is well defined. It is a scalar for spatial coordinate transformations. 
Notice that it is not purely local as it involves a quantity, J 72 , defined by the solution of an 
elliptic equation (|3.7|) . 

In practical calculations of the total gravitational energy, ()3.12j) is not so good because 
one needs to solve equation l|3.7|) . We now give an expression for Eq in terms of metric 
components and their derivatives. That expression does not depend on the existence of a 
solution of equation (J3.6|) . To this effect, we first notice that calculations near infinity are 
often simpler in spherical coordinates rather than x, y, z or axial coordinates R, z, <p like in 
section 4 (Hi). If the metric of the limiting flat space has a general form 

ds 2 = (dx ) 2 - j kl dx k dx l , (3.13) 

it is useful to take ds 2 as a flat "background metric". 14 This is done by replacing in IZ^i = 
lkmlirJZ mn partial derivatives by V^-covariant derivatives in the 7-background and the 

13 The original paper |10| contains precise statements about this theorem and useful references. It gives an 
example of an axially symmetric space in which the solution of (13.711 does not exist. The reason it does not 
exist is that TZ is not positive everywhere. 

14 The method that consists in introducing a second metric avoids the non-covariance and allows the use 
of coordinates that are not Minkowski coordinates at spatial infinity. The method goes back to a paper of 
Rosen I3ti|. Details on how suitable background metrics and mappings between the backgrounds and physical 
spacetimes may be introduced, in particular in asymptotically flat spacetimes, can be found in see also 

m 



13 



7-Christoffel symbols IT] (7) by the tensors 

A H = r Ti(l) ~^ki = \ 7 mn (Vklm + Vllnk ~ V„,7h) • (3.14) 
T's are Christoffel symbols in 7-space. Thus we write IZki as follows: 

K M = V m A£] - V k AZ + A£]A™ „ - A&A^. (3.15) 
Now remember a familiar identity usually used to construct Einstein's Lagrangian: 15 

U = C-V l k l where ^ = 7-^(77*') and £ = - 7 fc/ (A£}A™ „ - A&A^) . (3.16) 

We then substitute 1Z given by (|3.16j) into (|3.2|) . The asymptotic conditions on the metric 
have the result that, like in ()3.11|) . 



1 

2k 



Vik l dV = -Mc 2 . (3.17) 



Thus (|3.2|) may now be written as 

Eg = ~I {C + f[\B 2 + B-(Ax 8)]}dV. (3.18) 

ZK Jx° 2 

dV, as should be remembered, is not the proper volume element on x° = const [see (j3.1|) ], C is 
a scalar density because the A's are tensors. However, their value depends on the mapping on 
the flat background which is equivalent to choosing special coordinates. Eq can be calculated 
with this expression in terms of the metric components and their first order derivatives and is 
correctly given by the integral. Notice that if, for instance, Eq is found to be negative, then 
gravitational energy is negative irrespective of the choice of coordinates and irrespective of 
the nature of the sources of gravity. The formula is valuable in this respect. (|3.18|) is another 
form of the general formula obtained in [^5] (see the Appendix). 

It is at first surprising that the integrand does not contain derivatives of / because in the 
Newtonian approximation it is proportional to the square of d^f. But linearization is usually 
done in harmonic coordinates which, in particular, imposes (with Minkowski coordinates in 
the background): 

d l (^9 9 kl ) = 0, or equivalents -2^(^77^) = = 8 k . (3.19) 

This is how d^f reappears in Eq in the Newtonian limit. 



^Notice the first Vfe and the second Vs 
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It is also possible to write the total gravitational energy in terms of field components and 
a covariant integrand with the non-locality "shuffled into" a matter tensor integral. To that 
effect consider (|2,3U[) and (|2,31j) in the following combination: 

lu R kl - Zf^Roo = 2kT% =K + AVi£ l - A£ 2 - ± fB 2 . (3.20) 

Replace 1Z by C — V^fc', while moving the £ and B terms to the left hand side: 

4£ 2 + - fB 2 + 2kT% = C + Vi(A£ l - k l ). (3.21) 

The volume integral of the divergence on the right hand side is equal to zero because of 
the asymptotic conditions as follows from (|3.11jl and (|3.17|) . Thus, after dividing by 2k, the 
integral of C reduces to 



— f CdV = - [ {2£ 2 + -JB 2 )dV + [ T^dV. (3.22) 
2k J k J 4 J 



2f 

Replacing the C integral in H3.18() , we obtain a formula for Eg in which the non-local character 

pk 
k 



appears in the matter tensor part Tt - 



E G = -- I [2£ 2 + - f(B 2 + B ■ A x £)}dV - [ TtdV. (3.23) 
k J x o 2 J x o 

In particular, for perfect static fluids (A = B = 0) with local pressure p (T^ = —3p): 

E G = -^-J £ 2 dV + 3 J pdV. (3.24) 

In the weak field limit, the first term represents twice the gravitational energy while, thanks 
to the virial theorem, the second term is minus once the same. 

The important formulas of this subsection are: (1) the gravitational energy density as 
calculated by observers at rest given in (|3.12jl in terms of a non-local field (which may not 
always exist), (2) the total gravitational energy expressed in terms of local fields given in 
(|3.18|) . (3) the total gravitational energy calculated in terms of local fields with a real density, 
the non-locality being given in terms of the source of gravity, (|3.23j) . 

(ii) Gravitational energy for observers moving orthogonal to constant time slices 

The mechanical energy is now denoted Em because it is different from Em- With w^, see 
(|2.35|) . and the metric, see (j2,36|) . 

E M = I T>fw u ^/=g~di:„ = [ T 00 fdV where dV = y^d 3 x. (3.25) 

JxO J X 
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Therefore with T 00 / = ±G 00 f and (I2~i2l . the gravitational energy Eq is now: 

E G = Mc 2 - i- | q (ft + K 2 - K kl K k ^j dV . (3.26) 

We may again eliminate Mc 2 and obtain the energy density 16 in terms of a non local field. 
Thus, by analogy with (|3.3|) we set 

ill = ^ x lu and T k = -d k x, (3.27) 

and demand that the scalar curvature of the conformal metric 1Z* = 0. This defines x by an 
equation similar to (|3,6|) or ()3.7j) : 

£ = -4V fc ^ fe + 2.F 2 or V 2 e* /2 - -TZe^ 2 = 0. (3.28) 

8 

The asymptotic condition on x is the same as ([3.8)1 . If = and the dominant energy 
condition (jl.8|) is satisfied, then T 00 > and thus TZ which is proportional to T 00 + K kl K k i 
is positive. It then follows, see below ()3.8|) . that a solution for x always exists, whereas we 
have not proved this for static observers and, instead of ()3.26|) . we obtain, 

E G = -- [ {T 2 - \ K kl K kl )dV = I e G dV , (K = 0). (3.29) 
K J x o 2 J x o 

Note that the KK term has the opposite sign to the J- 2 term. The quantity is the 
gravitational energy density in 7-space, i.e., on x° = const. For practical calculations of the 
total gravitational energy, we may use a procedure similar to that used in (|3.13|) - ([3,18|) . and 
obtain Eg in terms of local fields: 

Eg = ~I {C-K kl K kl )dV in which £ = £(7^7) , (K = 0), (3.30) 

and also a formula similar to ()3.23() : 

E G = -- [ (£ 2 - K kl K kl )dV - [ (T k + W k T°)dV , (K = 0). (3.31) 

K Jx° J 

16 Natario [SSI wrote lESOJ, using in the form V k £ k = £ 2 + § fB 2 - / _1 k;(Too - § gooT) and called 

£ 2 + ~ the field energy density because of its analogy with electromagnetism. This formal definition is 
unrelated to our ea except in conformastatic spacetimes with metrics ds 2 — f(dx") 2 — f~ 1 dr 2 , see below, 
where the two "densities" are indeed equal to £ 2 . 
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4 Is gravitational energy negative? 

In Einstein's theory, gravity may be repulsive. 17 In Newton's theory, gravity is always at- 
tractive and gravitational energy is negative. In relativistic gravity we know, for instance see 
|25| . that even in the linearized theory of gravity, Eq is negative when the source of gravity 
is a perfect fluid but we do not know if it is true otherwise. Nevertheless, there are important 
instances in which Eq < and some of these are considered here. 

(i) Static spacetimes and static observers 

In static spacetimes A = and, see Q2.25j) . 1Z = 2kTq. The energy condition (J1.8|) is here 
^oo = /Tq > and implies 1Z > 0; in this case, as already shown by Cantor 0, a solution of 
Q3.7|) exists and 

E G = E G = — f CdV = — ( T 2 dV < 0. (4.1) 



K J K 

In spherical symmetry, the 3-space is conformally flat and the metric is of the form 

ds 2 = f{dx ) 2 -bdf 2 . (4.2) 



/ and b are functions of r. The case has been studied in some detail in where the 
connection with MTW's formula for gravitational energy has been given explicitly. Using C 
to calculate Eq we readily find that 

E G = -- [{(b-V^fdV. (4.3) 

Since outside matter, 

»=( 1+ £) 4 - <«> 

any spherically symmetric stellar model of isotropic radius tm has a total energy 

CM 2 / m \ 

E star = Mc 2 - E G {r > r M ) = Mc 2 + = Mc 2 1 + -— . (4.5) 

2r M V 2r A / / 

In the limiting case of a Schwarzschild black hole, ru = m/2 at the horizon and the energy 
of the hole itself is IMc 2 . This is the "quasi-local" energy attributed to the black hole by 
Brown and York [S]. 

(ii) Conformastationary spacetimes for static observers 

17 This follows from Raychaudhury's equation |35| . See Ellis and van Elst |14| for a covariant presentation 
(on the web). Repulsive gravity is discussed in various places see, for instance, 
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In conformastationary spacetimes |38| . 7^ is itself conformally flat but with no special 
symmetry. If we replace / by e 2 ^, see ()2,26j) . the metric takes the following form: 

ds 2 = e 2 ^ (dx° - A k dx k ) 2 - e- 2x Y,(dx k ) 2 . (4.6) 

k 

1Z itself has the form of (|3.6|) and Eq is given by (|3.12|) in which J- 2 is now purely local. 

The special case in which x = ^ an d thus Tk = £k nas been the object of much scrutiny 
in empty spacetimes 18 and in spacetimes with sources |24| . The gravitational energy for these 
metrics has been studied in some detail in |25| : the following additional remarks are of some 
interest. 19 If T h = £ k , (ETT21 can be written in the following form which is slightly more 
complicated but, as can be seen, it shows that Eq is manifestly negative: 

E G = ~~j f[\(B + Ax£) 2 + \{A ■ £f + r\l - \fA 2 )£ 2 ]dV < 0. (4.7) 

It is manifestly negative because, see (|2.7|) . 

g oo = r l (1 _ M 2 } = detidkO > Q (4g) 

Expression ()4.7|) gives a purely local gravitational energy density when expressed in spatial 
conformally flat coordinates. 

It is interesting to rewrite the last expression in a form which is similar to the energy 
density of the electromagnetic field. With 

A * = ljl/2 A j B* = VxA* and £* = ~ V/ 1 / 2 , (4.9) 
()4.7I) has this form 

E G = -- [ {B* 2 + 4/- 1 [(^4* • £*) 2 + (1 - A* 2 )£* 2 ]}dV < 0. (4.10) 

K J 

There is a striking resemblence between Eq and the energy of the electromagnetic field 
Eem- The electromagnetic energy momentum tensor is, in standard notation: 

C = ^ f T + ^" F %) where F l w = d ll A v -d v A lt , (4.11) 

and so 

Eem = J V=g~(Tt?) EM w v d^ = J {T$) EM dV. (4.12) 

18 See I38| for a review of the subject and references to original works. Solutions belong to a class of 
Einstein-Maxwell equations. With VkA k = 0, the metric is in harmonic coordinates, see 
19 Beware of the definition of B in |25| which is slightly different. 
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The electric field components are E^ = F^o = dkAo while those of the magnetic field are 
B k = rj kmn d m A n . In those terms Eem takes the following form in our coordinates 1(2.3(1 : 

Eem = ^J{B 2 + tr'KA* ■ Ef + (I - A* 2 )E 2 )]}dV > 0, (4.13) 

where B 2 = jkiB k B l etc... . Notice that in ((4. 13(1 .4* is the gravitational "vector potential" 
defined in (|4.9|) . not the electromagnetic vector potential A. 

(Hi) Axially symmetric spacetimes and zero angular momentum observers 

In axial symmetry, the metric (|2.36|) takes this special form: 

ds 2 = f(dx ) 2 -a ( dxK ) 2 ~ h ^ ~ udx ) 2 ; (4.14) 

K=l,2 

f,a,b and u = d(p/dx°, the angular velocity of the zero angular momentum observers as 
measured at infinity, are all functions of x K (K, L, M, ■ ■ ■ = 1,2). Under these conditions, 
the second fundamental form components are: 

K LM = K 33 = and K L3 = ± f-^bdiu. (4.15) 

Consequently, the gauge condition K = holds and T exists; using we find 20 that 

E G = -^J (T 2 - -f^bVw ■ Vuo)dV. (4.16) 

A particular case which gives a simple expression for Eq is when the 7-metric is conformally 
flat, i.e. b/R 2 = a, and 

ds 2 = f(dx ) 2 - a \dR 2 + dz 2 + R 2 (dip - uidx ) 2 ] . (4.17) 

This metric may be written in the following form: 

ds 2 = f A (dx° - A A Rdy) 2 - a \dR 2 + dz 2 + (1 + a' 1 f A A A )R 2 dip 2 ] , (4.18) 
f A = f-aR 2 u 2 and A A = -f^aRxv, (4.19) 

which makes it clear that it is not a conformastationary metric. To calculate the total 
gravitational energy we better use ((3.30(1 with C rather than T and a 'background' metric 

ds 2 = (dx ) 2 - dR 2 - dz 2 - R 2 dif 2 . (4.20) 

2 



20 Notice that Vw-Vw = a" 1 J2 ( dKLV 

K=l,2 
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The gravitational energy is 

E G = -— /(Vlogo- Vloga- - f^Wu- Vu)dV. (4.21) 
2k J 2 

The ZAMOs see less kinetic energy since they are dragged in the direction in which the source 

rotates and, therefore, one expects gravitational energy to be less negative. It is surely the 

case for slowly rotating sources, and presumably the case for fast rotating ones. 

One may use (|3.30j) to calculate the energy of an axially symmetric distribution of matter 

and, in the limit, of a Kerr black hole. An expression for quasilocal energy of a Kerr black 

hole, based on paper has been given by Martinez jHJ. There exists a great number of 

quasilocal expressions for energy which give different expressions for the energy of Kerr black 

holes gj. 

Appendix: Gravitational energy and conservation laws 

A different derivation of (|3.18j) . closely related to classical conservation laws in general 
relativity, was given in [25] and the result appears at first strikingly different from what is 
given here. The derivation of the original formula is very straightforward if one deals with 
stationary spacetimes. 

We start from the old Freud 21 [T^] complex i?\ antisymmetric in the two upper indices. 
In the notations 22 of [25] where a background metric is used we have 

V\ = -gx P D a (g p[fl g u]r7 )- 

A remarkable property of the tensor is this |34j : multiply by any vector field, say, h x ; the 
divergence of the resulting anti-symmetric tensor density, a conserved or divergenceless vector 
field density, is of the following form: 

Tjf is the matter tensor and t iL v is Einstein's energy "complex" in Rosen's (HEj covariantized 
form. If h x is replaced by the field of, for example, observers at rest [see (|2.1|) ]. the flux of 
x w x through a sphere at infinity, as shown by Freud, is 

/ t^ x w x dS^ = [ Vlt° l odS l = Mc 2 . 

J_ J r^oo 

21 A contemporary accessible version on the web can be found in |34|. 

22 The notations are those of the present paper with some additions: multiplication by v / — g of a tensor 
like g^ v is represented by a hat like g^ v = y '—gg^ LV and, for a change, = g^ v ' / 'y '—g while D stand for 
4-covariant derivatives in the ^-background spacetime. 
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Since the volume integral of T^w v is equal to the mechanical energy Em, it follows that 
gravitational energy is now given by the following expression: 

E G = + ^ x D v w x )dE u = J (f Q + t 0l Q £^dV. 

It takes some more calculations to show that, if the metric is of the form (|2.3|) . this integral 
is indeed (pTTHj) . 
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